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Abstract 

An axisymmetric panel code and a three- 
dimensional Navier-Stokes code (used as an inviscid 
Euler code) were verified for low speed, high angle-of- 
attack flow conditions. A three-dimensional Navier- 
Stokes code (used as an inviscid code), and an axisym- 
metric Navier-Stokes code (used as both viscous and 
inviscid code) were also assessed for high Mach number 
cruise conditions. The boundary layer calculations were 
made by using the results from the panel code or Euler 
calculation. The panel method can predict the internal 
surface pressure distributions very well if no shock 
exists. However, only Euler and Navier-Stokes calcula- 
tions can provide a good prediction of the surface static 
pressure distribution including the pressure rise across 
the shock. Because of the high CPU time required for a 
three-dimensional Navier-Stokes calculation, only the 
axisymmetric Navier-Stokes calculation was considered 
at cruise conditions. The use of suction and tangential 
blowing boundary layer control to eliminate the flow 
separation on the internal surface was demonstrated for 
low free stream Mach number and high angle-of-attack 
cases. The calculation also shows that transition from 
laminar flow to turbulent flow on the external cowl 
surface can be delayed by using suction boundary layer 
control at cruise flow conditions. The results were 
compared with experimental data where possible. 

Symbols 


Cf skin friction 
C p pressure coefficient 
c chord length 
w corrected weight flow 
x x coordinate 
a angle-of-attack 
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Introduction 

The progress in computational fluid dynamics in 
recent years makes it possible to use CFD codes as part 
of the design process. So far the panel method with a 
boundary layer calculation 1-4 has been used as a design 
tool for subsonic inlets. It is an efficient and economical 
way to obtain a good solution for internal subsonic flow 
without separated flow and shock formation; however, 
this method is limited to the internal flow of the 
inlets. 5 "® A three-dimensional Navier-Stokes solution, on 
the other hand, can also be obtained for a subsonic 
inlet;® however, the CPU time of 50 hr or more for one 
calculation on a Cray-YMP make it impractical for use 
as a design tool. For a high angle-of-attack case, a three- 
dimensional Euler/boundary layer approach can be used. 
The CPU time on a Cray-YMP for this approach is 
about 25 min which is reasonable for design applica- 
tions. An axisymmetric Navier-Stokes code can also be 
used for an ADP nacelle at cruise with less than 15 min 
of CPU time on a Cray-YMP. Because of the effect of 
downstream conditions on the upstream flow field in 
subsonic flow, the best results are obtained by combin- 
ing the inlet and the nozzle as an integrated geometry 
and performing the computation on both at the same 
time. 

In this paper, the axisymmetric panel/boundary 
layer method and a three-dimensional Euler /boundary 
layer approach were applied to a low speed, high angle- 
of-attack case, and the axisymmetric Euler/boundary 
layer approach and an axisymmetric Navier-Stokes 
calculation were applied to an ADP nacelle at cruise. 
The computational results are compared with experi- 
mental data. 

Experiments 

An ADP inlet/nacelle has been tested extensively 
with a 17-in. powered fan simulator in the NASA Lewis 
9- by 15-Foot Low Speed Wind Tunnel (LSWT) and 
8- by 6-Foot Supersonic Wind Tunnel (SWT) as part of 
a joint effort between NASA Lewis Research Center and 
Pratt & Whitney (P&W). 10,11 The experimental data 


provided an excellent database for computational 
analysis verification. The teat configuration which was 
called the conventional inlet was selected for the com- 
parison of analytical results with experimental data. A 
sketch of this configuration is shown in Fig. 1. The 
conventional inlet was the longest of a series of inlets 
designed by P&W and tested in the NASA wind tunnels. 

Numerical Methods 

The panel code starts with a geometry program 
which creates the control points for computation. The 
incompressible potential flow program is then used to 
calculate the basic solutions to the problem. 13 These 
solutions are used to provide a solution that satisfies the 
inlet operating conditions of freestream velocity, angle- 
of- attack, and inlet mass flow. Next, the incompressible 
flow is corrected for compressibility effects. 1 The com- 
pressible potential flow solution is then used as an input 
to the boundary layer program 4 which calculates the 
laminar, transition and turbulent boundary layer 
characteristics which can include the effects of suction 
and tangential blowing boundary layer control. There 
are three versions of the panel code; namely, two- 
dimensional, axisy m metric , and three-dimensional. 

The Navier-Stokes codes, PARC3D and PARC2D/ 
AXI, 12,13 were used to obtain the Euler solutions by 
turning off all viscosity terms (in the Euler mode). Both 
codes use the Beam and Warming approximate factor- 
isation algorithm. This algorithm is an implicit scheme 
which solves the set of equations produced by central- 
differencing the Navier-Stokes equations on a regular 
grid. Since these equations are formulated in the strong 
conservation form for a curvilinear set of coordinates, 
the algorithm is quite general. The artificial dissipation 
terms were added for improving stability and reducing 
the oscillations in the solution due to central differenc- 
ing. It includes a multi-blocked scheme so that it cm 
handle complicated geometries. In addition, boundary 
conditions can be easily specified by the user such that 
it can be called a user-friendly code. 

The surface flow information of the Euler solutions 
were used as an input to the TRACEON code 3 which 
traces a three-dimensional streamline on the nacelle 
surface. The flow information on this streamline is then 
used as an input to the boundary layer program 4 to 
obtain the flow characteristics such as skin friction, 
laminar/turbulent transition, etc. 

Computational Grids 

The two- and three-dimensional computational 
grids were generated with GRIDGEN2D. 14 The two- 


dimensional grid was rotated with respect to the x-axis 
to create an axisymmetric three-dimensional grid. 
Although this report is for an axisymmetric conventional 
ADP nacelle, the three-dimensional codes are available 
for a three-dimensional nacelle. 

The computational grids for the Euler calculation, 
are shown in Fig. 2. Figure 2(a) shows the two- 
dimensional grid with four blocks for the axisymmetric 
calculation. Only 10 200 grid points were needed for the 
axisymmetric Euler computation. For the axisymmetric 
Navier-Stokes calculation, only the C-grid around the 
nacelle was regrided for the viscous computation, and 
the other three blocks were kept the same as the Euler 
grid and computed in the Euler mode. By this zonal 
method, the total grid points, even for the Navier-Stokes 
calculation, was kept at 13 030. Figure 2(b) is the grid 
with five blocks for the high angle-of-attack computa- 
tion. A C-grid was wrapped around the nacelle for better 
resolution in the area of high curvature near the high- 
light (the leading edge of a nacelle). An H-grid was used 
for the other part of the computation. A total of 120 000 
grid points were needed for the three-dimensional Euler 
grid. In order to obtain a good prediction of external 
static pressures from the leading edge up to the trailing 
edge of the nacelle, two blocks were used for both grids 
to compute the exhaust plume from the nozzle. 

Results 

Low Mach number, high angle-of-attack operation 

Figure 3 shows the pressure coefficients on the sur- 
face for a free stream Mach number of 0.2, an angle-of- 
attack of 25 deg and w = 17.4 kg/sec. Although the 
pressure peaks in the computation al results were slightly 
different from experimental data, the overaR prediction 
for the high angle-of-attack using both panel code and 
PARC3D code in Euler mode was fairly good. 

Figure 4 shows the prediction of separated flow on 
the internal surface at a Mach number of 0.2 and an 
angle-of-attack of 25 deg by using the j>anel/boundary 
layer method and the three-dimensional Euler/boundary 
layer approach. The separation point was predicted at 
the point where the skin friction coefficient, Cp was 
equal to zero. Both methods predicted the identical 
separation point at x/c = 0.264. 

By using either suction boundary layer control or 
tangential blowing boundary layer control, the flow 
separation on the internal surface could be eliminated 
for the high angle-of-attack case as shown in Figs. 5(a) 
and (b). Herring’s boundary layer code 4 was used for the 
computation. A very small region of suction, applied 
from x/c = 0.258 to x/c = 0.269 with a bleed rate of 
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0.05 percent of the captured mass flow rate, could 
remove flow separation entirely as shown in Fig. 5(a). 
The internal flow separation also could be eliminated by 
tangential blowing at x/c = 0.242 with a slot height of 
0.5 mm and a blowing mass flow rate of 0.12 percent of 
the captured mass flow rate as shown in Fig. 5(b). 

Cruise performance 

The pressure coefficients on the internal and exter- 
nal surfaces of the nacelle are shown in Fig. 6 for a 
cruise Mach number of 0.85. The results include the 
axisymmetric Euler (PARC2D/AXI) calculation and 
three-dimensional Euler (PARC3D) calculation. Excel- 
lent agreement between the Euler solutions and ex- 
perimental data was obtained except near the shock 
location. The shock location is evident in the pressure 
distribution at an x/c of about 0.264. In this calculation, 
the PARC2D/AXI code captured the shock better than 
the PARC3D code. 

Figure 7 shows the comparison of the pressure coef- 
ficients between the Euler solution and the Navier- 
Stokes results at the cruise condition. The analytical 
result from the viscous Navier-Stokes calculation com- 
pared more favorably than the Euler calculation. The 
shock was captured by both Euler and Navier-Stokes 
calculations. In addition, the Navier-Stokes result 
eliminated the overshoot predicted by the Euler calcula- 
tion at about x/c = 0.3. 

The transition from laminar flow to turbulent flow 
at cruise conditions can be predicted by finding the 
critical point where the displacement thickness Reynolds 
number of the boundary layer first becomes greater than 
the critical Reynolds number. Herring’s boundary layer 
program 4 was used. Based on this method, the transition 
point for the ADP nacelle at cruise conditions is located 
at about x/c = 0.08, where x is the distance from the 
highlight and c is the length of the nacelle. The transi- 
tion point could be moved further downstream to 
x/c = 0.35 as shown in Fig. 8 by using suction boundary 
layer control, applied from x/c = 0.006 to x/c = 0.097 
on the external surface. For this case, the suction mass 
flow rate was 0.7 percent of the captured mass flow rate. 

Concluding Remarks 

An axisymmetric panel code and a three- 
dimensional Euler code were verified for low speed, high 
angle-of-attack flow conditions by comparing with 
experimental data. The numerical results from an 
axisymmetric Navier-Stokes code in both viscous and 
inviscid modes compared favorably with experimental 


data for high Mach number cruise conditions. The 
boundary layer calculations were made using the results 
from the panel code or Euler calculation. These calcula- 
tions demonstrated that both suction and blowing 
boundary conditions can be applied effectively to 
prevent internal flow separation in the inlet. The results 
also suggested that the suction boundary layer control 
is an efficient way to delay the transition point on the 
external cowl surface at cruise. These results provide 
confidence in using these approaches as part of the 
design process for a nacelle of an ADP propulsion 
system. 
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(a) 2-D grid for cruise computation. 
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Figure 3. — Pressure coefficient, Mach number = 0.2, a = 25 
w = 17.4 kg/sec. 



Figure 4.— Prediction of flow separation, Mach number = 0.2 
a = 25°, w = 1 7.4 kg/sec. 


(b) 3-D axlsymmetric grid for high angle-of-attack. 


Figure 2. — Computational grids. 
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(a) Suction boundary layer control 
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(b) Blowing boundary layer control. 

Figure 5. — Boundary layer control for a high angle-of-attack 
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Figure 6. — Pressure coefficient at cruise, Mach number = 0.85. 
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